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Abstract 



We explore the problem of time in quantum gravity in a point-particle analogue model 
of scale-invariant gravity. If quantized after reduction to true degrees of freedom, it leads 
to a time-independent Schrodinger equation. As with the Wheeler-DeWitt equation, time 
disappears, and a frozen formalism that gives a static wavefunction on the space of pos- 
sible shapes of the system is obtained. However, if one follows the Dirac procedure and 
quantizes by imposing constraints, the potential that ensures scale invariance gives rise to 
a conformal anomaly, and the scale invariance is broken. A behaviour closely analogous 
to renormalization-group (RG) flow results. The wavefunction acquires a dependence on 
the scale parameter of the RG flow. We interpret this as time evolution and obtain a novel 
solution of the problem of time in quantum gravity. We apply the general procedure to the 
three-body problem, showing how to fix a natural initial value condition, introducing the 
notion of complexity. We recover a time-dependent Schrodinger equation with a repulsive 
cosmological force in the 'late-time' physics and we analyse the role of the scale invariant 
Planck constant. We suggest that several mechanisms presented in this model could be 
exploited in more general contexts. 
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I. INTRODUCTION 

One important difficulty in canonical quantum gravity is the problem of time [1-3], which arises 
from the refoliation invariance of general relativity (GR). The Hamiltonian form of general relativ- 
ity (GR) has a quadratic constraint at each space point and, at the classical level, many-fingered-time 
evolution. Canonical quantization in the case of a spatially closed universe, leads, at least for- 
mally, to the Wheeler-DeWitt equation, which does not contain time. This is the frozen-formalism 
problem. As yet, there is no agreed way in which our manifest experience of a flow of time can 
be matched to structure in the quantum theory. In addition, there are severe technical problems 
in implementing the Hamiltonian constraints by quantum operators. 

In this paper, we shall not address these technical issues but will suggest that the frozen- 
formalism problem might be solved simultaneously with a further, less often noted problem - 
the failure of GR to be exactly scale invariant. This is best explained in terms of the GR initial- 
value problem, for which the only robust and effective method of solution is York's based on 
foliation of spacetime by spacelike hypersurfaces of constant mean (extrinsic) curvature (CMC) 
[4-6]. Represented this way, GR is a theory in which conformal three-geometries interact with a 
single global variable, the rate at which space expands. This global degree of freedom, a last ves- 
tige of Newton's absolute space, is why GR is not exactly scale invariant. The failure is mysterious. 

It is possible to construct an exactly scale-invariant theory of conformal three-geometries in- 
teracting with themselves alone - essentially GR without the global variable [7]. There also exists 
a closely analogous point-particle model of that theory [8]. Its key feature is replacement of the 
standard Newtonian potential by one homogeneous of degree —2, which ensures scale invari- 
ance. Newtonian gravitational forces are exactly recovered but with additional attractive forces 
that ensure scale invariance but are vanishingly small except on cosmological scales. 

The observation on which this paper is based is that Dirac-style quantization of the theory 
with this potential leads to a conformal anomaly. This breaks the scale invariance and leads to a 
renormalization-group (RG) flow governed by a dimensionsless scale parameter. We show that 
this leads simultaneously to 'time evolution' governed by the RG scale parameter. 1 We call this 
double emergence and put it forward as a possible resolution of two fundamental problems: the 
frozen formalism and the lack of exact scale invariance in gravity. 2 

1 The idea of considering RG as the origin of time in Shape Dynamics is due to S. Gryb. See also [9]. 

2 Its predictions can also be compared with the model explored in [10, 11] with Newtonian (degree —1) potential in 
which the single global degree of freedom that breaks scale invariance is traded classically for an internal time, after 
which the remaining dynamical degrees of freedom are all scale invariant. There are in fact two ways to achieve 
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In this paper we give a self-contained treatment of the classical particle model with the degree 
—2 potential of ref. [8] and its quantization with double emergence of the RG scale and an effective 
time through breaking of the scale invariance. We also treat the three-body problem, which gives 
us an opportunity to introduce what appears to be an interesting distinguished initial condition 
'at the origin of RG time'. We also explore a further novel feature of the quantization: it has 
dimensionless Planck constant. 

A word about concepts. In the Newtonian three-body problem all 3 x 3 particle coordinates are 
classical and quantum observables. But if the three particles are the entire universe, the position, 
orientation, and size of the triangle they form cease to have any observable meaning. There is also 
no external clock to measure time. Only the possible shapes of the triangle remain as observables. 
The collection of all possible shapes is shape space S. It is our fundamental concept. We insist that 
only functions defined on S are classical or quantum observables. Position, orientation, and size 
are all gauge. But we shall see that relative size, without becoming an observable, can play the role 
of an independent variable like time. 



II. THE ANALOGUE MODEL 



General relativity in the ADM formulation is a reparametrization-invariant theory, so its 
Hamiltonian vanishes on the solutions of the equations of motion. The particle model of [8] has 
a similar dynamical structure to GR but on a different phase space: R 6N , the space of the Carte- 
sian coordinates of N point particles r J = (ri J , T2 J , r^ J ), or, interchangeably, r J = (r x J , r y J , r z J ), 
together with their conjugate momenta p J = (pi J ,P2 J ,P3 J ) = (Px J ,Py J ,Pz J ), where J = 1, N. 
The reduced phase space is defined by the algebra of constraints 

N J J 

H = y E_lE_ + v = o , (i) 

^ 2mj ' w 

J=i J 

which implements reparametrization invariance, and 

N N 

P = J> J = 0, L = ^r J xp J = 0, (2) 

J=i J=i 

which implement translational and rotational invariance [12]. If V is a function of the N(N — l)/2 
separations | |r J — r 1 ] \ only, the algebra of constraints closes: 

{Li, Lj] = eijkLf- , {Li,Pj} = eijkPk i {PiiPj} = , (3) 



scale invariance, one (which we study here) stronger than the other (studied in [10, 11] ). 
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N av - av 

w p <> = E qp = - ^} = E = > ( 4 ) 

J=l « J=l * 

where we used Einstein's summation convention and {•, •} denotes the Poisson bracket 

{r{,pj} = S i:j 5 IJ . (5) 

This defines a theory of the evolution of relative configurations that can still differ in size. To imple- 
ment scale invariance and define a theory on shape space, we add the constraint 

N 

D = Y,* J -P J - (6) 

J=i 

Following [8] we call D the dilational momentum. It has the same dimensions as the angular mo- 
mentum L and measures expansion of the system just as L measures rotation. A scale-invariant 
theory clearly cannot allow physical expansion, for that requires an absolute scale. Vanishing 
of (6) is a necessary condition for scale invariance. But it is not sufficient: (6) imposes further 
constraints on V. Indeed, the Poisson algebra now also includes 

N ov 

{A Pi} = --Pi, {D,L j } = 0, {H,D} = -2H + ^rj— j+2V, (7) 

j=i i 

and it closes if X^jLi r j§^j + 2 V = 0. By Euler's theorem this requires V to be homogeneous of 
degree —2 in the coordinates rj. The potential with this property chosen in [8] is 

t r ^New T/ V-v mimj 2 1 A II / J\\2 fo\ 

V =^T> ^New = -2^ \\ Y I_ Y J\\ > P = Jp 2^ rnjrmWr -r || , (8) 

KJ 11 11 KJ=1 

where M is the total mass and MR 2 is the centre-of-mass moment of inertia. The classical me- 
chanics of this model reproduces Newtonian gravity to high accuracy in small enough subsystems 
but with a long-range attractive cosmological force that prevents the universe from expanding [8]. 
Indeed R is a constant of motion. 

In this paper we will study the quantum mechanics of this model. It turns out to be rich in 
interesting features. We first discuss an old issue in the quantization of constrained systems. 



III. QUANTIZATION: BEFORE OR AFTER REDUCTION? 



There are at least two ways to build a relational quantum theory. At first glance it seems most 
natural to get rid of the coordinate redundancies immediately and to quantize directly the true 
degrees of freedom. The other possibility is to impose the constraints quantum mechanically: this 
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is Dirac quantization [13]. The two procedures are in general inequivalent, so which is the correct 
one? 

The motivation for quantization after reduction is that the physical degrees of freedom live 
in the reduced phase space and therefore all the remaining ones are redundant (gauge) and un- 
physical. However, there are severe mathematical difficulties, generally insuperable, associated 
with working with true degrees of freedom, and for this reason alone the redundant degrees are 
usually retained. 

We have four reasons not related to mere computational convenience to do the same: 

1. There are several examples in physics of anomalies, in which a classical symmetry is not 
preserved by the quantum theory and the consequences have been confirmed experimentally. 
Quantization is by no means a unique procedure, as already operator ordering issues show, and 
a single classical theory usually corresponds to the classical limit of several quantum theories. 
Ultimately the observations tell us which quantization procedure provides the theory that fits 
them better. 

2. In gauge theories the interactions are local when expressed redundantly but most awk- 
wardly nonlocal in true degrees of freedom. 

3. The laws of transformation from one set of redundant coordinates to another are very sim- 
ple. In fact, Cartesian coordinates, like many other examples of redundant coordinates associated 
with an underlying Lie group, are crucial in suggesting the form of the kinetic energy. The point 
is that there is no unique prescription for choosing true degrees of freedom and, accordingly, no 
indication of how much kinetic energy should be associated with their rates of change. In con- 
trast, the Cartesian coordinates are distinguished and democratic: weighted with the respective 
particle masses, they provide an obvious measure of the kinetic energy that should go with each. 
The residual ambiguity that arises from time-dependent group transformations between Carte- 
sian frames is eliminated by means of constraints. Moreover, this approach selects a preferred 
form of these constraints, whereas in general one could get classically equivalent but quantum- 
mechanically different theories simply by multiplying the original constraints by arbitrary regular 
functions. This ambiguity is resolved if one works with group- theoretically defined redundant 
coordinates, which are distinguished and select the usual form of the Hamiltonian constraint. 

4. The reduction before quantization approach also runs into the problem of recovering a 
meaningful notion of evolution in reparametrization-invariant theories like GR. 
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IV. QUANTIZATION OF THE ANALOGUE MODEL 
A. Eliminating the momentum and angular momentum constraints 

It is easy to show that the two quantizations are equivalent for the momentum constraint P = 

0, which can be conveniently solved by moving to mass-weighted Jacobi coordinates {p K }, K = 

1, N — 1. These are a relational iV-body generalization of the reduced-mass coordinates used to 
decouple the centre-of-mass motion in the two-body problem. There are many choices of Jacobi 
coordinates for one and the same iV-body system. We refer to the Appendix IX A for details. 

Suppose we have chosen Jacobi coordinates {p J }. We denote by {tt j }, tt j = (vri J , tt2 J , ^3 J ), 
the momenta conjugate to {p J }. In the new coordinates, the moment of inertia is mapped into the 
barycentric moment of inertia [14]: 

N-l 

R"=^P J -P J - (9) 

j=i 

The new coordinates automatically implement the constraint P = 0, so we just need to consider 
the constraints L = and D = 0. The dilatational momentum constraint (6) becomes 

N-l 

D=^ P J -7T J = 0, (10) 

J=l 

and the Hamiltonian and angular momentum constraints are [3] 

W-l J J rr N 

#=£^ + ^ = 0, L = £p J x^ = 0. (11) 

j=i j=i 

The constraint L = can be solved by a gauge fixing (see Appendix IX B). 3 However, our main 
focus is scale invariance. From now on, therefore, we will keep the rotational gauge unfixed and 
solve the constraint L = by Dirac quantization. An exception will be the three-body problem, 
which is always planar when L = and has a simple algebraic phase-space reduction. 

B. Quantization of the Hamiltonian constraint 

The Poisson brackets between p J and tt j are 

{pf,nj} = 6 ij 8 IJ . (12) 

3 The 50(3) gauge symmetry of iV-body dynamics is probably responsible for a lot of interesting features of this 
model, in particular at the quantum level, as shown by the wealth of results obtained by Littlejohn and Reinsch [14], 
who took seriously the fibre-bundle structure of the configuration space. 
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We quantize these relations by promoting p\ and ttj to operators p\, nj that act on the Hilbert 
space = Sf 2 (M 3iV_3 ) and close the Heisenberg algebra 



' i //m <\j S IJ 



(13) 



A caveat: h S i is not the familiar Planck constant, but a dimensionless number (pj has the dimen- 
sions of a length, while is a length -1 ). We will see in Sec. VI how may be related to the usual 
dimensionful Planck constant. 

The quantum Hamiltonian constraint (11) is 

1 



Hip 



2 4A 3 ;v-3 + ^p )r-u. 



(14) 



where A3jv-3 is the Laplacian in 3iV — 3-dimensional Euclidean space. 

In polar coordinates {p\ , p\, p\, . . . , pf ~ 1 , p%- 1 , pg~ 1 ) -> {R, a), R G K+, a G S 3N ~ 4 , 



Hi/j(R,a) 



3iV-4 9 



+ 



where A S 3iv-4 is the spherical Laplacian, and we call V^ape = R ^New the shape potential [8, 12]. 4 
Its very interesting properties are investigated in the companion paper [11]. 

The reparametrization invariance of a relational theory like ours leads to a time-independent 
Schrodinger equation (15), the analogue in particle systems of the timeless Wheeler-DeWitt equa- 
tion in geometrodynamics. 



+ 2 Khape) 



^(R,a) = 0, (15) 



C. Separation of the problem 

The potential V = V^ ew /R allows us to separate the variables in Eq. (15) 5 

ip(R,q) = Yl ^( R ) M*) ■ 
n 

Substitution of this ansatz in (15) gives, for each n, the eigenvalue equation 

-f^shape = (~^si A S 3N-4 + 2 V^hape) ¥n = A„ <p n ] 



(16) 



(17) 



4 Saari [15] calls V^hape the configuration measure. 

5 We shall show that the functions tp n are countable, so we adopt this notation already. 
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this is a "scale-invariant equation". The operator H s h ap e, which acts only on the shape variables, 
is self-adjoint on the dense domain H 2 (S 3N ~ A ) (Sobolev space of «Sf 2 functions whose derivatives 
up to the second order are Jz? 2 ). By the Weyl and Rellich-Kato theorems [17] and the compactness 
of S 3N ~ A it can be proved it has a bounded from below and discrete spectrum. 

We still have to impose the constraints Lj if) = and D ip = 0. The former acts only on the 
scale-invariant wavefunctions because it commutes with both R and d/dR. Lj commutes also 
with i^shape/ an d therefore the remaining constraint Lj ip n = just acts as a selection rule on the 
scale-invariant eigenfunctions tp n . Applying this selection rule gives solutions defined on shape 
space, that is M? N /Similarity (M 3 ). 

For each shape eigenvalue A ra , we get a radial equation on «Sf 2 (M + ): 

which can be nicely rewritten by means of the change of variables u n (R) = (fiR)^ N -^l 2 £ n (R), 
where fx > is an arbitrary scale with dimensions length -1 , introduced to keep u n dimensionless: 

fl£d Un = + §«n = , 5n = A„ - 4 (3iV _ 4 )(3iV - 6) ^ . (19) 

D. The anomaly 

The anomalous behaviour of a 1/i? 2 potential like the one in the radial equations Eq. (19) has 
been well-known since the pioneering work of Case [18] in the 1950s. Initially dismissed as non- 
physical, its importance has only relatively recently been recognized in very different areas, from 
molecular physics to black holes (see e.g. [19] and the references therein). In particular there has 
been a surge of interest after the experimental confirmation [20, 21] of the Efimov effect [22], which 
is now better understood as a consequence of this quantum anomaly. 

The novelty of our approach is that we found, starting from relational first principles, similar 
equations describing a cosmological setting. A difference is that each radial equation is coupled 
with scale-invariant degrees of freedom and ultimately, through the anomaly, will determine their 
evolution in an emergent time. Let us develop the argument. 

6 One of the authors has proposed, in [16], that a notion of time might be hidden in the frozen wavefunction resulting 
from the reduction before quantization approach. This might be the case if the wavefunction of the universe is 
concentrated on 'time capsules', i.e., configurations which suggest a history. In the reduction-before-quantization 
approach, the only allowed solution of (17) is one with X n — 0, corresponding to zero total energy and a frozen 
wavefunction. Rewriting (17) in the form flshape ¥>n = —A S 3N-A<p n + (2/7i|)K,hape¥'n = (An/fiJ) ip n , we see that 
plays the role of a dimensionless coupling constant. It might perhaps be possible to adjust its value to ensure the 
existence of the zero eigenvalue. This would be a novel first-principles derivation of the strength of quantum effects. 
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We can see that a 1/R 2 potential can generate an anomaly in this way: consider the quantiza- 
tion of the dilatation constraint (6). There is a mild ordering ambiguity. Without loss of generality 
we will take the Weyl ordering 

/> = _£Wtf + tfA0 = -_(i«^ + j ^). ™ 

J,i v 7 

From Eq. (7) and the canonical quantization we get 

/ JV-l Qy\ 

\ J=l "Pi J 

Classically the quantity A vanishes for a 1/R 2 potential, due to Euler's theorem, but quantum 
mechanically this theorem breaks down, as we will see. The anomaly can be rewritten: 

R 3N-5 / pV 



A = 2F + p-V(F-) = 2^-V-^^j , (22) 
where p = (p{, . . . , p{ , . . . , p^" 1 ) and V = (d/dp\, d/dp{, d/dp^ 1 ). 

A ~ 2 ^~ y " \rW^1 J ~ 2~ ' ) ^hape + J^P ' W shape . 

Defining ^3at_4 as the surface of a 3N — 4 dimensional sphere, we have 

V • (^3) = V m ^5 m - 3 (p), A = R 3N - 5 n 3N . 4 5 3N - 3 (p)V shape + r^p ■ VV sl . (23) 

Equation (23) is exactly where the classical Euler theorem departs from the quantum one. Now 
£3iv-3^ _ _____ ^- R3N _ 4 S(R), and because V_hape does not depend on R Euler's theorem holds for 

it. This implies 

A = V shape R- l 5(R) . (24) 

The distribution A is called an anomaly. Our algebra does not close even if the potential is homo- 
geneous of degree —2, 

[H, D] = -2H + ihV^R-^iR) , (25) 

unless we ask for particular boundary conditions at the origin. We can interpret this result as a 
failure of Euler's theorem at the quantum level [19, 23]. In the next section we will see that the 
closure of the algebra is incompatible with self-adjointness of the radial Hamiltonian. We should 
emphasize that self-adjointness is a working hypothesis that we retain in a theory of the whole 
universe, for which it is not a priori clear whether all the basic postulates of quantum theory, 
established in subsystems of the universe, are preserved. Caveat lector. 
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E. Self-adjoint extensions 



The radial equations (19) on the dense domain 



S>(H^ ad ) = {u, u absolutely continuous and u(0) 



u'(0) = 0} 



(26) 



have been studied extensively in the literature [19]. The deficiency indices depend on the energy 
regime. They are 

1. (0,0) if 5n/^i > 3/4, so that the radial Hamiltonian is self-adjoint on Jz^ 2 (IR+). 

2. (1, 1) if —1/4 < <7n/^si < 3/4, so that there is a U{\) group of self-adjoint extensions, 
labelled by the index 9 . Scale invariance is generally broken, but there are two conformal 
fixed points of the renormalization group (RG) flow. We will call this the weak coupling 
regime. 

3. (1, 1) if 5n/^si < — l/4, but each element of the U(l) group of self-adjoint extensions 
breaks scale invariance, i.e., there are no conformal fixed points in the RG flow. A limit 
cycle behaviour emerges. We will call this the strong coupling regime. 

From the fundamental characterization of self-adjoint extensions [17], the domain of Hq is the 
direct sum of @{H™ ad ) and the span of 



consisting of a linear combination of the solutions of H™ ad u = ±i [iu; we need fi for dimensional 
consistency; 1Z = fiR/h S i, and E n is the order of the Bessel function Kg n . The self-adjoint exten- 
sions break scale invariance, generating 'bound states' [19] with nonzero eigenvalues. With H@ 
denoting the ^-extension of the radial operator H™ ad , the eigenstates are given by the equation 



9n /h 2 si + 1/4 , 



(27) 




d 2 u g n u k 2 




(28) 



&R 2 h 2 si K 2 " fi 2 



The substitution u(JZ) = \Jn1Z/ nf(7l) leads us to a modified Bessel equation for /; the two inde 
pendent solutions are by definition 




(29) 



but, for both real and imaginary values of H n , only the first one is a solution in J£ 2 (M). 



12 



To choose a self-adjoint extension amounts to imposing a boundary condition at the origin in 
the previous problem, so we now want to understand which eigenfunctions sjnll/fj, K-= n (jjTljr 
corresponding to eigenvalues — k 2 , are compatible with the boundary condition 

lim u(K) = lim w e (K) . (30) 

This condition gives a relation between k and 9. Expanding we and u near the origin, we get 

„-2 



exp 



1 fcos(0/2-^~ n /4) 
— log 



(31) 



> cos(0/2 + 7rE n /4), 

The weak and strong coupling regimes have qualitatively different behaviours. In the weak 
coupling regime E n G [0, 1] is real and Eq. (31) is one-to-one, so each self-adjoint extension selects 
one and only one state of the discrete spectrum. The function n n (9) is defined only in the domain 
9 G [0, 7r — 7r H n /2] U( 7r + 7T S n /2, 27r], where it is bijective (see Fig. 1, left-hand side). The two 
boundary points, ir — tt E n /2 and it + ir E n /2, corresponding respectively to the values k = and 
k = oo, merge at the transition S n — > between the weak and strong coupling regimes (Fig. 1, 
right-hand side) . This is relevant for the RG interpretation of the problem, as we will see later. 

In the strong attractive regime S n is imaginary, and this causes an abrupt change. The function 
k(9) is defined on the whole range 9 G [0, 2ir), but is multivalued: in this case Eq. (31) can be 
rewritten as 



exp 



+ 



arctan 



sin I 



2-7T m 



m G N . 



(32) 



cos 9 + e 71 "! 5 "!/ 2 , 

Each 9 selects an entire discrete spectrum made of a tower of states with an accumulation point 
at zero. These are all the values of the energy compatible with the boundary condition fixed by 



K/fl 




3„ = 0.6 



A 



2/r 



s»= o 



*-9 



2/r 



FIG. 1: k vs 9 in the weak regime. The two preferred choices of 9 correspond to k = and n = oo. 
The boundary between the weak and strong regimes, where the two distinguished self-adjoint extensions 
merge, is at H„ = 0. 
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FIG. 2: k/h vs 9 in a cylindrical representation in the strong regime. At each 9 there is a conformal tower 
of states in a geometric sequence: the ratio between two neighbouring levels is a constant fixed by 3 n . 
Continuous variation of 9 over a 27r range, leads from one level of the tower to the next. 

the chosen 9. The motion of each level as 9 is allowed to vary leads to a limit cycle behaviour, as 
explained in Fig. 2. 

If we imagine the line twisting around the cylinder as a spring, a change in H n corresponds 
to a change in its elongation. In fact the ratio between one level and the one immediately 
below, is equal to a constant that depends on E n : E^/E^^ = e 71 "/^™'. This is known as 

the Efimov effect [22]; it is characterized by the presence of a 'conformal tower' of states in a 
geometric sequence {..., (^) e - 27r /|Sn| ; (« a .) e -5r/[ s "[ ) (^eVI^ (« a ) e 27r/[H„[ ) j 

V. RENORMALIZATION GROUP TIME 
A. The Callan-Symanzik equation 

The anomaly broke the constraint Dip = 0, which does not need to be satisfied by the wave- 
function anymore. But there is a domino effect involving the constraint H ip = 0, which turns out 
to be incompatible with the requirement of self-adjointness of the radial Hamiltonian. 7 We will 



7 This is just a consequence of the fact that the eigenvalue zero in the strong regime is an accumulation point in the 
spectrum of every self-adjoint extension of each radial Hamiltonian without being an eigenvalue associated with a 
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shortly argue that this generates an RG flow. 

The minimal ansatz that makes the theory consistent is to relax the Hamiltonian constraint by 
introducing an energy eigenvalue — k 2 , 

HiP = -k 2 ?P, (33) 

which is not in conflict with the other constraints P ip = and L ip = 0. Support for this mini- 
mal ansatz is that it maintains the reparametrization invariance of the classical theory though no 
longer with energy eigenvalue E = 0. Moreover from Eq. (31) this requirement fixes one and 
only one self-adjoint extension for each radial equation (19) at the subtraction point. They are the 
extensions 6 n such that 

n nm {dn) = ~k 2 Vn . 8 (34) 

In other words, the single ansatz (33) not only makes the theory self-consistent but also fixes all 
the mathematical redundancies that arise from the theory of self-adjoint extensions. 

Let us now study how the notion of time emerges in our model as an RG flow. The wave- 
function ip splits into a sum in which the coefficients u n will be interpreted as couplings for the 
scale-invariant eigenfunctions <p n , 

1>(R, a) = u n {R)R^ N ^ 2 vn{<r) , a G 5 37V " 4 ; (35) 

n 

the couplings u n flow differently under the RG flow, as detailed in the following subsection, in a 
way that is controlled solely by the scale-invariant energies g n with which each is associated. In 
other words, each scale-invariant eigenfunction determines its own RG evolution. Moreover, from 
(35) we see that an evolution of the u n 's amounts to a redistribution of the relative amplitudes 
of the scale-invariant eigenfunctions. This redistribution looks very similar to an evolution with 
respect to a 'time' that is the scale -parameter of the RGflow. 9 

This might appear rather formal, but the language of effective field theories provides us with 
a physical interpretation of what is going on. The Hamiltonian constraint implies Eq. (28) for the 
u n 's. This equation can be recast in an interesting form if we re-express it as an equation for the 

proper eigenfunction. 

8 The label m does not affect the choice of 8 n in the strong regime because it corresponds to a complete turn around 
the cylinder (2). Moreover in the weak regime the label rn is not present at all. 

9 An advantage of this proposal is that the usual internal time scenarios generically have great difficulty in meeting 
the requirement of monotonicity whereas in our framework time, being associated to RG flow, is monotonic by 
definition. 
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logarithmic derivative of u n : 



r> du n (R) 

ln(R) = ^gr- , u n (R) oc ef'° sR ^(R')dio S R> (36) 
u n (R) 

The radial equation for j n becomes 

Uln 2 , i 9n i r>2 2 /or,-, 

^ „ = -7n + 7n + Z2" + R K ■ ( 37 ) 

a log i? rtgj 
This last equation has the form of a Callan-Symanzik equation for a coupling ^ n with beta- 
function 

/?(7n, fl) = = - (in ~ 1/2 - v^ + W) ( 7n - 1/2 + ^ + i2 2 « 2 ) , (38) 

Kaplan a/. [24] and Mueller et al. [25, 26] have already appreciated the equivalence of the 
Schrodinger equation for the 1/R 2 potential and the RG flow equations for a theory with this 
particular beta function. Following [25, 26], we provide an interpretation of the running of the 
couplings: Eq. (28) has a bad singularity at R = (a total collision, the 'big bang' of the model), 
and therefore one will seek to regularize the potential near the origin. 10 For example one could 
take a a square-well cutoff [19]: 

{q n I R? , R > Rq 
(39) 
1 /Rl , R < R 

Then one has to match the solution and its derivative at the cutoff, and an efficient way to do this 
is to match the logarithmic derivatives at R = Ro: 

lim — 22- = lim — ?g- . (40) 

r ^r+ i>(R) r^r- i>{R) 

This equation implies a relation between the coupling 7 and the cutoff radius Ro, which causes 
7 to run. If we calculate the logarithmic derivative of the wavefunction inside the cutoff radius 
at Rq, we obtain 7 n (i?o)- As a function of Rq, 7 satisfies Eq. (38) [24]. The radial Schrodinger 
equation is equivalent to the RG group flow equations (37). 

The usual interpretation in RG theory is that changing the cutoff corresponds to probing the 
system at different energies or, equivalently, at different scales. Actually the independent variable 
in an RG flow is not an energy but a ratio of energies. This matches the intuition that, if we are 

10 Regularizing the potential is a complementary approach to that of self-adjoint extensions because, like it, it amounts 
to fixing a boundary condition near the origin. But the self-adjoint extensions approach is non perturbative and 
regulator-independent, so we keep both languages at hand to supplement each other and guard ourselves against 
artefacts of the regularization choice. 
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considering the entire universe, there is no scale external to it. In our cosmological model the 
RG flow parameter is a ratio of scales, and it gives rise to a running of the couplings u n 's that 
redistributes the probabilities among different scale-invariant eigenf unctions. The Universe at 
different (relative) scales will change, as the RG flow guides the probability from one region of 
shape space to another. 11 



Let us sketch the evolution of the couplings u n under the RG flow. This determines the relative 
weights of the wavef unctions in (35). There are three kinds of behaviour: 

1. None of the u n 's corresponding to g n > 3/4 evolve. They are frozen at their initial values. 

2. All the u n 's corresponding to — 1/4 < g n < 3/4 evolve under the RG flow between two fixed 
points of (38). One of them, which fixes k = 0, is a Friedrichs extension that effectively sets 
the anomaly to be zero. The other one corresponds to n = +oo. The u n 's have the following 



Because the u n 's depend only on the product k(8)R it is clear that there is a duality between 
the evolution under the RG flow and the evolution of R, from the definition k(8)R = R{9). 
The radius of the universe, measured as a dimensionless ratio relative to the initial (sub- 
traction) point, changes under the RG flow, renormalizing the couplings u n : 



The two fixed points are reached for R = and R = oo. The infinite past and future 
are conformal: this is a realization of Strominger's proposal [27] that the evolution of the 
Universe results from the renormalization of a conformal field theory flowing between two 
fixed points, associated to dark-energy dominated epochs of accelerated expansion. 12 The 
structure of this model, however, is conceptually richer: we also have a frozen regime and 
a cosmological Efimov effect. 

We will not discuss here the major interpretational issue - which as yet has not definitive answer - that every quan- 
tum theory of thel Universe must face: what is the meaning of these probabilities? 
See also the model proposed recently in [28] and related papers. 



B. Renormalization-group flows 



expressions: 




(41) 




(42) 
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3. If g n < — 1/4 we enter the Efimov effect region. The most general expression for each u n is 
a linear combination of all the eigenfunctions of the extended Hamiltonian: 

u n (R(9)) = V^m(0)nK En (K m (e)TV). (43) 

m 

We have the same correspondence between evolution under RG flow and evolution of R 
as in the weak regime. Moreover, as we saw, the limit cycle behaviour implies that the 
following geometric relation holds: n m = e m7r / = ™ k. 



VI. SEMICLASSICAL LIMIT 



In this section we will study how the physics of our toy model appears at late times, far from 
R = 0. In this limit we can regard R as a "heavy" degree of freedom, in a way that we will make 
precise in the following discussion. 13 Without loss of generality, let us rewrite the Hamiltonian 
constraint (33), substituting ift(R, a) = r](R)x{R, cr). The evolution under RG flow is not unitary, 
so the function i] is fixed by the requirement that the function x stays normalized on shape space, 
(X, x) = 1/ where the scalar product is the shape Hilbert space Jz? 2 (£ 3Ar_1 ). Denoting with a prime 
the derivatives w.r.t. R, 

~ ftsi(v"x + 2 VX + X"v) - H si 3N R 4 ( r l'x + VX) + Jp{-ftlAs3 N ~±X + ^ s hapeX) = -^VX ■ (44) 

Denoting by (-, •) the scalar product on the scale-invariant Hilbert space Jz? 2 (5 3JV-4 ), and assum- 
ing that x is normalized to 1 at each value of R, so that {x,x} = 1/ the mean over the scale-invariant 
("light") degrees of freedom of (44) gives 

-ny - 2/&/( x , x') - hjrtx, x") - 4^F^ W + v(x, x'» + (45) 

Jj2 i- k si(X, A 5 3iv-4X) + (X, ^shapeX)] = -^V (X, X) = ■ 

Let us suppose now that the mean values over the light degrees are vanishingly small, apart from 
(^shape) = (X; ^shapeX)- This is very reasonable because of our assumption that the shape degrees 
of freedom are very rapidly varing w.r.t. R. We get 

- ^si 7 /' - R V + "IpP" 7 ? = ■ (4 6 ) 

13 We do not have the "multiple choice problem" of usual semiclassical approaches because any function / of R and 
D, f(R, D), is singled out by the anomaly. This is a significant advantage of the shape-dynamic approach. 
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Now substitute again this equation in (44) to get 

- 2h 2 si HL\ X ' - h 2 si X " - ^si 3 ^ 4 ^ + ^(-h 2 si A s3N . iX + KhapeX - (^shape)x) = . (47) 

Under the hypothesis that 77 is a complex field, it can always be rewritten as t]{R) = 
a(R)e lS ( R ^ h si. 14: We substitute this in (46) and require the real and imaginary parts to vanish 
separately 

fc2 // , c'2 t 3iV — 4 , (Vshape) 2 n fAQ\ 

— n si a + a b — fi S i — a H — ^ — a + k a = U (48) 

R R 

-2h si -S'- h si s" - h si S' = (49) 
a R 

Let us examine at Eq. (48) in the classical limit defined compatibly with our previous assumptions 
as the case in which the amplitude a varies much less than the phase S. In this limit Eq. (48) is a 
Hamilton-Jacobi equation 

y 2 = _ K 2_ (5Q) 

R £ 

From S' = dR/dt em , the Hamilton-Jacobi equation gives us a definition of an emergent classical 
time; it has the dimensions of a surface, 

JO 

df m = . _ (51) 

-. (^shape^ 
~ l Wn2T- 

Let us recast physics in this emergent classical time. Substituting (50) and (49) in (47) 

S" 1 

- 2 i k si S' X + ^si ^7 X ~ ^si X" + ^2 ( _/i si A 53^-4 + V^hape ~ (^shape)) X = (52) 

At late times we can disgregard the term S" /S' oc 1/R?. The term in x" measures the kinetic 
energy associated to the overall expansion of the Universe. If this is neglibible w.r.t. the kinetic 
energy associated to the all whole of the shape degrees of freedom, we get the equation 



ih, 



d X 1 



S1 dt em 



^2 (^si A 53iv-4X + (Khape ~ (Kshape))*)- (53) 



Because R is a heavy degree of freedom, -^A S 3n-4X ~ A R 3iv-3X/ the previous equation finally 
becomes a time-dependent Schrodinger equation: 

■k d X t2 a , 1 (^shape) 

tksi dt^ = ~ h n A ^ N - :iX + ^^NewX p2 ( 54 ) 

This equation is very interesting because it shows what could be achieved in a more sophisticated 
implementation of the ideas of this toy model: 



The assumption that 77 is complex could be unjustified because is taken to be a solution of a real equation (see [29]) 
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1. It proves that our time-dependent Schrodinger equation could arise (for late-time physics, 
in a semiclassical limit) from a completely scale-invariant cosmological model as a conse- 
quence of the scale anomaly. A cosmological force appears. It is proportional to the mean 
value of the scale-invariant shape potential and repulsive. 

2. If we ignore the cosmic origin of the dimensionful time (51) and consider (54) as the equa- 
tion of an effective theory empirically determined, we would regard time as a fundamental 
unit and, for dimensional consistency, h$i would be given the dimensions of an energy times 
a time. This is why we would write (54) with H instead of h s \, where h is considered a dimen- 
sionful constant. This suggests a general mechanism by which dimensionless fundamental 
constants in a scale-invariant model (in this case, h s \) can be replaced by dimensionful con- 
stants in a semiclassical approximation. 

3. The fact that the 'gravitational constant' depends on the 'time' R is a serious potential de- 
fect of the model that we do not wish to hide, especially since there are relatively strong 
bounds on the secular variation of Newton's constant. It should be remarked, however, 
that we derived Eq. (54) under the assumption that R is almost constant compared with the 
observable degrees of freedom, so we cannot use it for predicting a cosmological evolution 
of G. Moreover, we are mainly interested in this toy model as a conceptual test ground for 
exploring the problem of time and scale invariance. 

VII. THE THREE-BODY MODEL 

An important advantage of the three-body case is that the non-anomalous constraints P = 
and L = can be explicitly reduced before quantizing. But the reason why we are studying this 
problem is for the conceptual insights we will get from it. First we will get a physical under- 
standing of what the scale-invariant Planck constant means. Secondly we will propose a natural 
solution to the question of setting initial-value conditions in our scale-invariant model. Both is- 
sues are interesting because seems they will not be restricted just to our particle toy-model. 

A. Phase-space reduction 

We will follow Montgomery [30] in the phase space reduction. The first simplification comes 
from the fact that, because the angular momentum L is zero, the problem is planar, so we can 
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gauge fix two out of three angular momentum constraints plus one translational constraint by 
fixing the plane in which the motion stays. We are left with three two-component vectors r , 
r 2 and r 3 that define the position of the three bodies on this plane. Then we gauge fix the two 
remaining translation constraints by going to mass-weighted Jacobi coordinates: 



1 I m 1 m 2 , 2 i\ 2 m 3 (m 1 +m 2 ) ( 3 m 1 r 1 +m 2 r 2 \ 

p = \ — r — ( r - r ) , p = \ — z — z — r z • ( 55 ) 

y mi + 1712 V THl + m-2 + m 3 \ TXl\ + 7712 / 

We are left with four coordinates p 1 , p 2 and a single angular momentum constraint to gauge fix. 
The coordinates 

^^(H^H'-IIP 2 !! 2 ) . ^ = P 1 -P 2 , w 3 = p'Ap 2 (56) 

are invariant under the remaining rotational symmetry and therefore give a complete coordinate 
system on the reduced configuration space. The norm of the vector w is proportional to the 
moment of inertia 1 1 w 1 1 2 = R 2 /4, so the angular coordinates in the three-space (wi,W2,w 3 ) coor- 
dinatize shape space, which has the topology of a sphere [30]. This shape sphere, as can be seen in 
Fig. 3, has several distinct regions of interest. 

B. The spectrum of the scale-invariant Hamiltonian 

For our purposes it suffices to say that the scale-invariant Hamiltonian takes the form [30] 15 

-^shape = ~^si ^S 2 + ^shape , (57) 

where A S 2 is the 2-dimensional spherical Laplacian and 

V sha e = - V {munjfl 2 1^ ^ 

SKaPe i<j ( mi + m j) 1/2 ( m ~t +1^2 + m 3 ) 2 yj\ - cos x cos(V> - ipu) ' 

X is the elevation angle from the equator and tp the azimuthal angle on the shape sphere, ipu is 
the azimuthal angle of the two-body collision between particles / and J. Moreover, our Hilbert 
space will be Jz? 2 ^ 2 ), with inner product 

(f,g)= dsmxdi}f(x,^)g(x,^)- (59) 
Js 2 



15 Strictly speaking, the scale-invariant Hamiltonian is not self-adjoint due to the three singularities at the two-body 
collisions, and a U (3) family of self -adjoint extensions would be needed. But these singularities are not as bad as 
the 1/R 2 singularity at the origin, because they go like 1/R. Moreover, these singularities reproduce the Newtonian 
gravitational interactions among particles, and there are good reasons just to take the Friedrichs extension for all of 
them, which are the extensions that allow us to best model the hydrogen atom. 
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FIG. 3: The shape sphere of the three-body problem. The equator corresponds to the collinear configura- 
tions; at three points on it there are two-body collisions, where V^hape is singular. Between them, at positions 
that depend on the mass ratios of the particles, there are saddle points of T4hape (Euler configurations). The 
points in the southern hemisphere represent triangles which are mirror images of the corresponding ones 
(with the same longitude and opposite latitude) on the northern hemisphere. The north and south poles 
are equilateral triangles, at which, by a result due to Lagrange, K^ape attains its absolute magnitude for 
arbitrary values of the masses. The critical points of V^hape (Euler points and the equilateral triangles) are 
central configurations and play an important role in three-body theory. For example, if the system is 'held 
at rest' at one of them and then 'released', it will fall homothetically - without changing its shape - to a 
triple collision at the centre of mass, beyond which the solution cannot be continued. The figure is for the 
equal-mass case, with the binary collisions evenly spaced at 120° intervals. 



We will provide a qualitative description of the spectrum of i? s hape- We already mentioned 
that it can be proved that the spectrum is discrete and bounded from below. Let us start from 
the most excited states, with positive scale-invariant energy -f/ghape^n = <Pn and A n S> 1. A 
convenient complete orthonormal basis in which to calculate the matrix elements of -ff s hape i s that 
of the eigenfunctions of the spherical Laplacian, the spherical harmonics: 



-A s2 Y e m = £(£ + l)Yr. 
It can be shown that in the expansion 

(Y e m \H shape \Yf) = 5 mm ,5 ul l{l + 1) + / dsinxdV Y e m Y™' V si 



(60) 



(61) 



the matrix elements of V^ape are bounded from above, being V S i bounded and harmonic functions 
normalized, while the term £(£ + 1) grows arbitrarily large with £. Therefore asymptotically the 
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spectrum is well approximated by that of the free particle on the sphere, that is 

A n -> frgi ^ + !) > for n ' ^ sufficiently large. (62) 

For the less excited states, one gets a hyperfine splitting from the matrix elements (^ m |14hapel^ n ')/ 
which introduce an m-dependence of the eigenvalues. 

Another regime in which we can say something qualitative about the spectrum is close to the 
ground state. In fact our problem is completely equivalent to that of a particle constrained on 
a sphere and subject to a potential that has three Coulombian wells. Close to one of the three 
singularities \ = 0, tp = ipu, the potential is well approximated by a Coulomb potential and the 
Laplacian is approximately flat: 
const 

V shape -====== , A S 2 ~ + —2 , x, i> in a neighborhood of (0, ipu) . (63) 

\/X W-iplJ) dx "V 

If fr S i is sufficiently small, the ground state and the first few excited states will resemble those of the 
hydrogen atom (or, rather, a linear combination of those of three independent hydrogen atoms). 
The reason why h S i needs to be small enough is that h s i plays, for those wavefunctions, the role of 
a Bohr radius, setting a scale of quantum phenomena in shape space. Therefore it measures how 
far the wavefunction ventures from one of the singularities at (0, ipu). If it is localized enough, it 
will feel neither the curvature of the shape sphere nor the presence of the other singularities. 

In the limit of small K s i we are able to calculate the ground scale-invariant energy: it is (let us 
assume, for simplicity, that the three particles have the same mass) 

Ao~--^. (64) 
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Then from Eq. (19) we get 



90 ~ - Q + ^(3iV - 4)(3iV - 6)) ^ = - i- , (65) 

which can be made smaller than — 1/4 if h S i is small enough. This means that some shape eigen- 
states enter the strong coupling regime. 

C. The problem of the initial conditions in a scale-invariant model 

We have, at this point, a qualitative picture of the features of our model. At least in the three- 
body case, we have established the form of the spectrum of the scale-invariant Hamiltonian, esti- 
mated the ground-state eigenvalue and determined its asymptotic behaviour for large eigenval- 
ues. This spectrum is almost all that one needs to determine the i?-evolution of the wavefunction. 
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A further bit of information is needed, in the form of an overall normalization constant of each 
u n function. In fact choosing the self-adjoint extension (33) for all n only fixes the (logarithmic) 
derivative of u n in Eq. (28), which we reproduce here for convenience, 

~ h l ^W + W Un = ~ K2un ' 9n = \n-\m- 4) (3JV - 6) h 2 si . 

This is a second-order differential equation and therefore depends on two initial conditions. The 
logarithmic derivative v! ju is fixed Vra by the self-adjoint extension via Eq. (30), and the other can 
be fixed by specifying an initial condition for the functions u n , which can be set at a finite radius 

u n (R ) = c n , c n £<C , (66) 

or at Rq = (through a limiting procedure), or even at Rq = oo. 16 

We propose now a way of fixing this initial condition; it relies on the highly asymmetric struc- 
ture of shape space. In fact, C = — V^hape/ which may be called the complexity, defines a positive- 
definite function on shape space that takes its minimum value on the configuration of the system 
that is more uniform than any other configuration. In the case of the three- and four-body prob- 
lems for arbitrary values of the masses [31], the minimum of C is at the equilateral triangle and 
regular tetrahedron respectively. Batty e et al have done numerical calculations for up to 10 4 parti- 
cles which showed that the minimum of C is realized on remarkably uniform (super-Poissonian) 
configurations [31]. 

The fact that all relative configuration spaces have a distinguished point was noted by one of 
us many years ago [16, 32, 33] but before the importance of quotienting wrt dilatations as well as 
translations and rotations was fully appreciated. The distinguished point identified in [16, 32, 33] 
corresponds to the configuration with all particles collapsed to a single point and, for obvious 
reasons, was called Alpha. The additional quotienting wrt dilatations and the associated passage 
to shape space makes it possible to transform the qualitative intuition of the original idea into 
a quantitative theory because C is an objective scale-invariant measure of complexity on shape 
space. It is a classical and quantum-mechanical observable. Simultaneously - and suggestively - it 
is, as we have shown, (minus) the potential that governs the scale-invariant dynamics. 

We now define Alpha to be the most uniform shape of the system, where C has its minimum. 17 

16 We noted earlier, in footnote 6, that in the case of quantization after reduction the spectrum, if it exists at all, consists 
solely of a zero eigenvalue. If h si can be 'tweaked' to ensure existence of such an eigenvalue, the theory will still be 
incomplete if the eigenvalue is degenerate and corresponds to a superposition of zero-eigenvalue wavefunctions. 

17 There is no corresponding notion of a most complex state: there is an Alpha, but no Omega [8, 11, 34]. 
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The presence of the Alpha point in shape space suggests the existence of a distinguished initial 
condition for the wavefunction of the universe: at early times it could be peaked around Alpha. 
When this choice is made, the expectation value of K-hape is close to its maximum, the complexity 
close to a minimum. 

VIII. DISCUSSION AND CONCLUSIONS 

In this paper we have discussed the so-called problem of time in quantum gravity using an 
analogue particle model whose algebra of constraints is similar to the one of GR in Hamiltonian 
form. This model is almost uniquely fixed by relational first principles identical to those that one 
of us proved to be at the basis of GR [12, 35]. 

As we mentioned in the introduction, our basic idea is that two problems in GR could end 
up providing the solution to each other. We explored this possibility in the particle model. The 
two problems are the frozen nature of the Wheeler-DeWitt equation and the failure of GR in its 
usual representation 18 to be fully relational: a single global degree of freedom, the volume of the 
Universe, cannot be considered as gauge. In this paper we insisted on strict classical relational 
principles by choosing a fully scale-invariant analogue model. We then discovered that time and 
scale emerge naturally through an anomaly that arises if we attempt a consistent Dirac canonical 
quantization. 

We concluded that in Dirac quantization of the analogue model global scale invariance is an 
anomalous symmetry that is broken in the quantum theory. This is due to the singularity that the 
potential shows at a total collision. Its renormalization leads to an RG flow in which the radial 
eigenfunctions u n in the expansion 

tP(R, a) = ^ u n {R)R^ N ^/ 2 ^(a) , a € S 3 ^ , (67) 

n 

behave like RG running coupling constants for the shape eigenfunctions ip n . In particular their 
logarithmic derivative satisfies a Callan-Symanzik equation. The behaviour of the couplings u n 
in Eq. (37) depends on the scale-invariant energy g n associated to ip n (see Eq. (17) and (19) for the 
definition of g n ). Above a certain limit, the couplings are frozen; then there is a regime in which 
they flow between two conformal fixed points; and finally the lowest energies exhibit an Efimov 
behaviour with limit cycles. 

18 The caveat "in its usual representation" relates to the possibility noted in footnote 2 of achieving scale invariance by 
trading a global degree of freedom for an internal time. 
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In this intriguing picture the dimensionless ratio of the scale of the universe relative to a refer- 
ence value becomes 'time' in the sense that the changing scale of the universe renomalizes the cou- 
plings u n , and this corresponds to a flux of probability on shape space. The scale ratio works this 
way because its change produces an RG flow that redistributes the relative probabilities among 
the shape eigenfunctions. We call this proposal double emergence, because scale and time emerge 
from the same quantum phenomenon, an anomalous symmetry breaking. 19 

What we have produced is a concrete implementation of the holographic cosmology scenario 
(albeit for a particle toy-model) that has very sound principles attached to it. In fact, our motiva- 
tions have nothing to do with holography at all. This could be interesting from the point of view 
of finding a physical foundation of the AdS/CFT correspondence. 

Besides this main result, in this paper we obtained a number of other interesting results and 
observations, which we summarize here. 

1. The semiclassical regime. We studied the wavefunction of the system at large values of our 
'time' R, under the assumption that the wavefunction fluctuates much more rapidly in the 
shape directions than in R (Born-Oppenheimer approximation). With a WKB ansatz, we 
were able to find an effective time-dependent Schrodinger equation satisfied by the wave- 
function, with a function of R and (V^, a pe) playing the role of time. This effective semi- 
classical dynamics is expected to become dominant at late times (large i?'s). Our equation 
contains a repulsive cosmological force and a dynamically fixed Newton constant G = 1/R, 
but the equation is not valid if R varies at cosmological scales. 

2. The dimensionless Planck constant. The 'Planck constant' that appears in the canonical 
commutation relations of our model is a dimensionless number. This is due to the fact that 
the action is dimensionless, and the momenta have dimensions of inverse lengths. In the 
semiclassical, 'late-time' regime we were able to show a mechanism in which the familiar, 
dimensionful Planck constant could emerge. At the fundamental level time has no mean- 
ing, and the quantity that emerges at the classical level as time has the dimensions of a 
length squared. If one attributes to it a fundamental role, as is done when dealing with the 
effective physics that is relevant in the laboratory, then the Planck constant takes on some 
dimensions, which are needed to convert from squared lengths to times. This role of the 
Planck constant is identical to that of Boltzmann's constant in thermodynamics: fundamen- 



Admittedly this model is still affected by the quantum-mechanical measurement problem. The wavefunction in fact 
spreads over macroscopically distinguishable configurations in shape space. 
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tally it is dimensionless, but at the effective level one attributes to it the dimensions that are 
necessary to convert from energy to heat, which at the fundamental (in that case, micro- 
scopical) level are the same thing. Even if eventually the particular toy model considered in 
this paper gets ruled out, we expect such a role for the Planck constant to emerge whenever 
one assumes that the fundamental degrees of freedom are dimensionless. 

3. The three-body problem. We studied the case of three bodies, which is unique for the 
important simplifications that it allows: the angular momentum constraint can be alge- 
braically solved, and shape space can be conveniently described as a two-dimensional 
sphere, whose radius is proportional to the moment of inertia R, our 'time'. The shape 
Hamiltonian is simply the spherical Laplacian on the shape sphere plus V^hape/ which is a 
relatively simple function on the sphere. A qualitative study of its spectrum was enough 
to establish that all of the three regimes of the RG flow are realized. Moreover the scale 
invariant Planck constant appeared as a scale of quantum phenomena on shape space. 

4. The initial condition. The only thing that this model leaves unspecified is the initial condi- 
tion for the wavefunction on shape space, which can be fixed by setting the coefficients u n 
at a certain value of R. Developing the classical proposal of a distinguished 'a-point' in the 
configuration space of the scale-invariant degrees of freedom made in [34], we suggest that 
the driving potential of the scale-invariant dynamics, Vghape/ is also a natural measure of 
complexity in the system. In fact its maximum is attained the most uniform state, which we 
called the 'a-state' . 20 So a natural initial condition for the wavefunction on shape space is 
such that all the probability is concentrated on this uniform shape, the equilateral triangle 
in the three-body problem. The evolution would then make the wavefunction spread out 
of that point, driving change in shape. 

An interesting possibility is to extend these ideas to more general contexts. In particular it is 
intriguing to ask if GR could be considered as a fully conformal gauge theory on conformal su- 
perspace (the geometrodynamic counterpart of our shape space) in which a single global gauge 
degree of freedom becomes dynamical and a preferred notion of cosmological time emerges. The 
theory of shape dynamics [35-39] could be the natural context in which these ideas can be ex- 
plored. 



See VII C for a discussion which generalize this proposal to N body problems. 
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IX. APPENDICES 
A. Jacobi coordinates 

A simple algorithm [40] for creating Jacobi coordinates requires specification of a complete 
directed binary tree, whose leaves represent the N particles, and the N — 1 internal vertices are 
associated with the Jacobi coordinates. The algorithm works this way: given such an arbitrary 
tree, one assigns to each "parent" vertex four numbers: the three coordinates of the centre of mass 
of the two "child" vertices and the sum of their masses, which are taken of course from the leaves. 




FIG. 4: The binary tree algorithm for the three-body system: the two Jacobi coordinates associated with the 
internal nodes are p 1 = r 1 — r 2 , p 2 = (miY 1 + m 2 r 2 )/(rai + m 2 ) — (mir 1 + rn 2 r 2 + m 3 r 3 )/(mi + m 2 + m^). 

Then one defines the Jacobi coordinates as a function on the internal nodes. This function 
associates to each node the difference of the coordinates of its two "child" nodes: pP arent = r left — 
r right . When applied to the root node, this function gives the iV-th Jacobi coordinates, that is, 
the coordinates of the centre of mass of the system, which has now been decoupled and can be 
discarded, so that one can work with just with the first 3N — 3 coordinates. The great advantage 
of the coordinates, besides the centre-of-mass decoupling, is that they leave the kinetic metric 
diagonal [40]. Its eigenvalues in this coordinate system are the Jacobi effective masses fij, which 
are just the masses of the nodes associated with each Jacobi coordinate. This algorithm represents 
a relational way to move to the centre-of-mass reference frame. 

B. Diagonalization of the inertia tensor 

We can solve the angular momentum constraint with a gauge-fixing, that is, with a choice of 
orientation of our axes. A global gauge-fixing, though, has to satisfy certain requirements to be 
good. The main requirement is that the constraint surface defined by the gauge fixing must never 
be parallel to the vector field generated by the constraint it is supposed to gauge fix. This trans- 
lates into the requirement that the Poisson bracket between the gauge fixing and the constraint 
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must be invertible, which means it must be everywhere (weakly) non-vanishing. 

A popular choice of axes among iV-body specialists is the one defined by the principal axes of 
the moment-of-inertia tensor [14]. It is attractive for its simplicity and symmetry and more than 
anything else for its relational character, but it fails to be global for the reason mentioned above: 
it becomes degenerate at certain points of configuration space. 

The moment-of-inertia tensor is defined as 

N N 

jab = ^ mi {r[ - r c ™){rl - rg m ) - 5 ab ^ mj | |r 7 - r cm | | 2 . (68) 
i=i i=i 

Its eigenvectors are the three principal axes of inertia of the iV-body configuration. The matrix 
I ab is symmetric, and therefore can be diagonalized with a rotation that puts the three princi- 
pal axes of inertia respectively on the x, y and z axes. The three constraints that impose such a 
condition are the ones that set the off-diagonal elements to zero: 

/ 12 = j23 _ jl3 _ q ^ (69) 

This condition, though, does not always uniquely fix the axes for two reasons: 

1. When the system is in a symmetric configuration (spherically or axisymmetric) there are 
identical eigenvalues (two if axisymmetric, three if spherically symmetric), and the matrix 
is already (block) diagonal. 

2. Even in the non-symmetric case, the solution to equation (69) is not unique: there are three 
solutions, corresponding to the residual symmetry under exchange of the axes. But this is 
harmless as one can resolve the ambiguity by requiring the eigenvalues to be arranged in 
order of magnitude: I 11 < I 22 < I 33 . 
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